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We study travelling fronts of equations of the form utt + 4>(u)u x = u xx + f(u). A criterion for the 
transition from linear to nonlinear marginal stability is established for positive functions <f>(u) and 
for any reaction term f(u) for which the usual parabolic reaction diffusion equation ut = u xx + f[u) 
admits a front. As an application, we treat reaction diffusion systems with transport memory. 
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PACS numbers: 



I. INTRODUCTION 



Reaction diffusion equations are used to model trans- 
port phenomena in a variety of contexts such as popu- 
lation dynamics, transmission lines, flame propagation 
among others. The prototype of such equations, for 
which a thorough understanding (IJ ^ of its proper- 
ties has been achieved is the parabolic reaction diffusion 
equation 



ut 



u xx + f(u), with /(0) = /(l) = 0, (1) 



where subscripts denote derivatives. For positive reac- 
tion terms, sufficiently localized initial conditions evolve 
into a monotonic decaying travelling front joining the sta- 
ble u = 1 to the unstable u = equilibrium points. For 
bistable reaction terms, which satisfy / < in (0, a) and 
/ > in (a, 1), with J Q f > 0, it is possible to find initial 
conditions for which the system will evolve into a mono- 
tonic decaying front joining the two stable equilibrium 
points u = 1 to u = 0. In the first case there is a con- 
tinuum of speeds for which there exist monotonic fronts, 
the system evolves into the front of minimal speed. In 
the second case there is a unique speed. 

The use of Eq. (|l|) to model a physical process involves 
assumptions on the stochastic process that describes the 
motion of the individuals. More specifically, Brownian 
motion is assumed ||. If instead a more realistic pro- 
cess is considered then, in one spatial dimension, the 
differential equation that describes the motion is a hy- 
perbolic reaction diffusion equation of the form 



u tt + (f>( u ) u t = u xx + f(u). 



(2) 



ft has been shown that, as in the usual reaction diffu- 
sion equation, for positive 4>(u), the hyperbolic equation 
admits monotonic decaying fronts of speed Ch < 1 and 
that the speed of the front is determined from a related 
parabolic equation 

Equations of this form have been studied for particu- 
lar cases of <j>(u) obtained assuming specific functions for 
memory effects in the diffusion term |j| 0, ||, [§. Recent 
work has dealt with the stability of the traveling front 



JlOj l . For constant 4>(u) and positive, concave, / the 
traveling fronts are stable to local perturbations. The 
asymptotic behavior of the front of minimal speed has 
also been found . Numerical solutions for positive but 
non concave / indicate that the transition from linear to 
nonlinear marginal stability occurs at the same param- 
eter values than for the parabolic equation Eq. 
Detailed analysis for a piecewise linear reaction term has 
given insight into the nature of possible travelling waves 
0. In systems with exponential transport memory 4>(u) 
is of the form 



<j>(u) 



/'(«)• 



For this case a variational principle for the front of min- 
imal speed was constructed for positive reaction terms 
/(it); it was proved that for concave reaction terms, as 
in the parabolic equation, linear marginal stability holds 

In this work we construct a variational principle for 
any positive function 4>{u) and general reaction terms for 
which the parabolic equation admits a front. From this 
principle upper and lower bounds are constructed which 
permit the a priori determination of the transition from 
linear to nonlinear marginal stability pi [l^] . We recover 
all known results as particular examples, and construct 
others that show different types of behavior depending on 
the explicit functions f(u) and <f)(u). Our main result is 
the following, consider Eq. (2), with <f)(u) > 0, and f(u) 
a reaction term for which fronts of the parabolic equation 
exist, as described above. The minimal (or unique, for 
the bistable case) speed of the fronts joining the stable 
to the unstable points u — 1 to u = is given by 



c 2 



sup 2a 



tig(K(u))f(u)du 

So -g 2 (u)/g'(u)du 



where 



and 



a 



— = 4>{u)du 



K{u) = a (f)(u)du. 



(3) 



(4) 



(5) 



2 



The supremum is taken over all positive decaying func- 
tions g(u) for which the integrals exist. From this result 
it follows that, for positive reaction terms which satisfy 
/'(0) > 0, the following bound holds, 



> 2 (0) 



< 



1 



< 4 sup 



cj){u)K(u) 



which enables one to characterize functions for which 
marginal stability holds. The lower bound, the marginal 
stability value, is that imposed by linear considerations 
alone. Improved lower bounds valid for all reaction terms 
are obtained by direct use of the variational principle. 



II. SPEED OF FRONTS 

Consider the hyperbolic reaction diffusion equation 

u tt +(j>{u)u t =u xx + f(u), (6) 

with 



(u) > 



and 



/(0) = /(l) = 0, 



where f{u) and 4>{u) € C 1 [0, 1]. We shall assume that / 
belongs to the class for which monotonic fronts joining 
the equilibrium points u = 1 to u — exist. The precise 
conditions have been spelled above in the introduction, 
and in more precisely elsewhere ^ |l3j . 

We wish to find the minimal or unique speed for which 
there is a monotonic decaying travelling wave solution 
u = u(x — c*t) ofEq.(||). The speed satisfies (l — c 2 )u zz + 
c(j>(u)u z + f(u) = 0, limw^-oo = 1, limw^oc = 0, 
u z < 0, where z = x — ct. It is known that when / > 
there is a continuum of fronts for a range of speeds 
c* < c < 1, where c* is the minimal speed of a related 
parabolic equation Q. There is a unique speed in the 
bistable case. We may, as usual, consider the trajectory 
in phase space by defining p = —u z (u). The monotonic 
decaying front obeys 

(l_ c 2) p%--af>(u)p + f(u)=0, (7a) 
du 

p(0)=p(l) = 0, p>0 in (0,1). (7b) 

Before going any further we recall the constraints posed 
by linearization around the fixed points. Linearizing 
around u — we find that p approaches zero as mu, 
where m is the positive root of 



(1 - c 2 )m 2 - ccf)(0)m + /'(0) = 0, 



that is, 



1 



W(0) 

2 1 - c 2 2(1 - c 2 ) 



1 -c 2 



When /'(0) > 0, since it is known that monotonic fronts 
exist with c 2 < 1, m is real if the term in square brackets 
is positive, that is, 



1 



> 



4/(0) 
2 (O) 



(9) 



where we have called cl the minimal speed derived from 
linear theory or, as it is known, linear marginal stability 
value. 

When /'(0) < 0, linear theory imposes no constraint 
on the speed. Analysis near the equilibrium point u = 1 
imposes no constarint on the speed as we assume that / 
is positive or of the bistable type. 

The simplest, but not unique, method to obtain a vari- 
ational characterization of the speed is to introduce a 
stretching of coordinates which reduces the equation to 
the standard parabolic reaction diffusion equation ||. 
Since <p(u) > 0, we may introduce a new independent 
coordinate defined by the following transformation, 



V = K{u) 



(10) 



with K(u) defined in Eq.(^|). The transformation is in- 
vertible, and y varies between and 1. In the new coor- 
dinates, Eq.(R) reads 



dp 

p(y)-r - 



dy a(l — c 2 ) 



p(y) 



f(K-\y)) 



a(l - c^(K-^y)) 



= 0, 



(11) 



p(0)=p(l) = 0, p>0 in (0,1). 



This equation is the equation for fronts of the parabolic 
equation, of speed 



and a reaction term 



F(y) 



o(l - c 2 ) 



f{K-\y)) 



We now check that F(y) satisfies all the requirements of 
the existence of fronts, if / does. Effectively, F(0) = 
F(l) = 0; for <j> > and 1 - c 2 > 0, sign(F) = sign(f), 
and 



F{y)dy 



1 



1 -c 2 



f(u)du. 



So that if / is of the bistable type, for which fronts exist, 
the same holds for F if 1 — c 2 > which we know to hold. 
Moreover 



F'(y = 0) 



1 



a 2 (l-c 2 ) 



/'(0). 



3 



We may then apply directly the variational principle for 
fronts of the parabolic reaction diffusion equation |jl3), 
that is, 



c = max 



2 Jo 9(y)F(y)dy 
3 Io(-9 2 (y)/9'(y))dy 



(12) 



where the maximum is taken over all positive decaying 
functions g in (0,1) for which the integrals exist. The 
maximum is attained for 

9 = 9 = exp ^- J ~^ ds ^j > < yo < 1 (13) 

provided that F'(0) < 0, or, if F'(0) > 0, when m/c > 
1/2. This is, 

m 1 



. >-[l + ^c?-4F'(0)]. 
c 2 

Written in terms of the original quantities, c, /, <f>, this 
condition reads 

c 2 > 4/'(0) 



1 - c? (j) 2 (0) ■ 

As in the usual parabolic case the maximum is attained 
except at the linear marginal stability value. Here too 
one can show that the linear marginal stability value is 
obtained by taking the supremum instead of the maxi- 
mum. In terms of the original functions / and (f> , Eq. fll2] ) 
reads, 



1-c 2 



2af o \g(y)f(K-\y))/0(K-\y)))dy 

f \-9 2 (y)/g'(y))dy 



(14) 

Notice that the simpler form Eq. (^) is obtained by chang- 
ing to the original independent variable u. 



A. Upper bound 

We present here a simple derivation of the upper 
bound, valid both for positive and bistable reaction 
terms. The following evident inequality holds, 



Wo 1 



1 g(y)f(K- 1 (y)) , 



1-c 2 



LltlA -I 

9 f (-9 2 (y)/9'(y))dy 



< max sup 

9 y 



2af(K- 1 (y)) 



y4>(K-i(y)) 



S l y9{y)dy 
Jo -g 2 (y)/g'(y)dy 



But, as shown in the Appendix, 



l ,i 

yg(y)dy < 2 / -g 2 (y)/g'(y)dy, 
o Jo 



from where we obtain 

,,2 



C 



< sup 



II 



4a/(K" 1 (y)) 



y<t>(K-\y)) 



(15) 



Going back to the original independent variable, this is 

r 4a/(u) 



1-c 5 



< sup 



K{u)<j)(u) 



(16) 



Observe that this upper bound coincides with the linear 
value when the supremum occurs at u = 0. 



B. Lower Bounds 



III. CRITERION FOR LINEAR AND 
NONLINEAR MARGINAL STABILITY 

For the parabolic reaction diffusion equation Eq(|l|) , for 
/ > 0, we know that the minimal speed, cpf of the 
propagating front is bounded below and above. When 
these two bounds coincide, then one can determine the 
speed unambiguosly. A first estimate comes from the 
bound | 



2VT(0) < c PF < 2 sup 



/(«) 



When these two bounds coincide linear marginal stability 
holds. It may still hold if these two bounds do not coin- 
cide. This can be decided making use of the integral vari- 
ational principle which improves the lower bound 
and of a minimax variational principle which improves 
the upper bound 1 14 . Combined use of the two permits 
the exact determination of the speed. For the bistable 
case, we can obtain the speed making use of the integral 
variational principle. In this section the analogue of the 
above results is obtained for the hyperbolic fronts. 



The lower bounds have already been obtained, one 
is the linear bound Eq.(Q), and improved lower bounds 
can be obtained by choosing particular trial functions g. 
Equality, in the expression below, holds for a certain g. 



> max 



0n r 1 g{y)f(K- 1 {y)) rh: 
2a Jo 4>(k-Hv)) y 



g 2 (y)/g'(y))dy 



(17) 



IV. APPLICATIONS 

A. Constant <j>(u) 

As a first example consider the case (f>(u) = 4>o, a con- 
stant. As explained in the introduction the stability of 
the traveling front for this equation has been proven re- 
cently |[o]]. Numerical investigations have shown that 
the transition from linear to nonlinear marginal stability 
occur at the same parameter values as for the parabolic 
equation with the same reaction term @] . This last fact 
follows directly from the bounds above. For a constant 



4 



(j>(u), we obtain a = l/(f>o and y = u. Then, the bounds 
read 



4/'(0)<-^<4sup 

1 C u 



/(«) 



The upper and lower bounds are identical to those for 
the speed of fronts of the parabolic equation parabolic 
equation; the speed is given by 



2 , ,2 



1-C 2 



2 Jo 1 g{u)f{u)du 

tiax — 

5 ti(-9 2 (u)/g>(u))du 



C 



PFi 



(18) 



The speed of the front is determined in terms of the 
speed of the parabolic front with the same reaction term, 
therefore the transition from linear to nonlinear marginal 
stability occurs at the same parameter values as for the 
parabolic case. The speed itself, however, is lower than 
the speed of the front of the parabolic equation, 



i> + Cp F 



< c 



PF- 



(19) 



B. Suppresion of nonlinear marginal stability 

As a second example, which shows that substantially 
different behavior may occur, take the reaction function 
f(u) — u(l — u)(l + au). The front of minimal speed for 
the parabolic equation is given by 

C PF = 2 for < a < 2 
^2 in 



CpF 



for a > 2. 



If we now consider the hyperbolic equation with (f>(u) = 
1 + au, the upper and lower bounds coincide and the 
speed 

for all a. 



l-c 2 
Effectively we obtain 



Z 1 

C 1 — U 

4 = 4/'(0) < - < 4sup- - = 4. 



1 -c? 



1 + au/2 



The minimal speed is the linear marginal stability value, 
no transition from linear to nonlinear marginal stability 
occurs. A significant slowdown of the speed takes place. 



C. Time delayed diffusion 

In systems with time delayed diffusion the function 
(u) adopts the form 



form which arises when considering diffusion with expo- 
nential time delay. Here s is a parameter. For this func- 
tion a = 1/s and K(u) = u — f(u)/s. The speed of the 
front was studied for positive reaction terms, and it was 
shown that for concave reaction terms marginal stability 
holds H . This can be proved easily from the present ap- 
proach. As Eq.(|l6|) shows, linear marginal stability holds 
for all positive reaction terms with /'(0) > provided 
that the supremum of 



B(u) = 



4a/Qu) 
K{u)(j){u) 



(20) 



occurs on u — 0, since B(0) = 4/'(0)/c/f> 2 (0), is equal to 
the lower bound. To prove then that linear marginal 
stability holds it suffices to show that B'(u) < for 
u G (0,1]. Taking the derivative we obtain, after some 
algebra, 



B'(u) 



K{u) 



ff" sa 
4> 2 (u) 4>{u)K{u 



■h(u) 



where hiu) — uf — f. It follows immediately that for 
concave reaction terms B'(u) < 0. For concave terms 
/" < 0, so the first term in the square brackets is nega- 
tive. And h(u) is negative for concave functions as well, 
h(0) — 0, and h' — uf" < 0. So we have recovered the 
result that for concave functions with f'(0) > 0, linear 
marginal stability holds. As it occurs for the parabolic 
case, it may hold for noncave functions as well, fact which 
can be determined by use of the variational principle. 

Consider again the reation term /(u) = u(l— u)(l+au) 
as in subsection B. In the parabolic equation transi- 
tion from linear to nonlinear marginal stability occurs 
at a = 2. Choose here s = 2 which guarantees 4> > for 
< a < (5 + v / 2T)/2 w 4.79 . We show in Fig. 1 that 
the transition from linear to nonlinear marginal stability 
occurs at a lower value of a than that for the parabolic 
case. We chose as a trial function g(u) = (1 — u) 1 /u 1 - 8 . 
The dot-dash line shows the lower bound imposed by the 
linear marginal stability value, the solid line the lower 
bound from the variational principle. The transition oc- 
curs at at a as low as 1.6 at least. The precise transition 
value can be found by using a better trial function. The 
dashed line shows the upper bound, which guarantees 
that linear marginal stability holds for < a < 1, values 
for which / is concave. 



D. A bistable reaction term 

The variational principle and the bounds obtained 
above hold for positive and bistable reaction systems, 
no assumption on the sign of / is made in the deriva- 
tions. We may therefore apply them to a bistable reac- 
tion term. We take as an example a case for the speed of 
the parabolic front can be determined exactly, 



4>{u) = s - /'(«) 



f(u) = u(l - u)(u — a) for < a < 1/2 



■5 



The unique value of the speed for which a front exists is 

cp F = —= - aV2 

If <f>(u) = 4>o a constant, the speed of the front is given 
by Eq([l9]), and the exact speed is known with cpf given 
above. 

For this bistable reaction term, with <j>(u) = 1 — f'(u), 
which is positive for < a < 1/2, with a simple trial 
function g(u) = 1 — it, all integrals are elementary, we 
obtain, using Eq.([L4h, 



a 18 - 49a + 14a 2 



- 49a + 14a 2 ' 

Numerical determination of the speed can be made with 
any desired accuracy by use of improved trial functions. 

We wish to point that, in recent work, a variational 
principle valid only for positive reaction terms has been 
applied to a bistable example ||. Bistable reaction 
terms, by definition, are not positive throughout the 
whole interval, so the hypothesis of their derivation is 
violated. We attribute the nearly perfect agreement be- 
tween the speed of the front obtained from their (non- 
applicable) variational principle to that obtained from 
direct numerical integration to a coincidence for the spe- 
cific trial function chosen. It is not difficult to find an 
acceptable trial function, sharply peaked at the origin, 
for which the integrand of their variational expression, 
and so the speed, becomes complex. 



which marginal stability holds. We have seen that de- 
pending on the nonlinearities of the equation the range 
of validity of linear marginal stability, can be increased 
or decreased, compared to the comparable range for a 
parabolic equation with the same reaction term. The 
stability of the fronts has been proved for a limited class 
of such equations [[To[ , numerical solutions have been ex- 
plored for systems with memory effect. These results will 
be more relevant if the stability of fronts of more general 
systems is proved. At present numerical results indicate 
that the front of minimal or unique speed can be selected. 
The general initial value problem for hyperbolic reaction 
diffusion equations is an open problem. 

APPENDIX A 



Let g{y) be a positive decaying function, which satisfies 
g(l) = 0, and call h(y) = —g'(y), so that h > 0. Then 



l -,2 

yg(y)du 



;(-l)(^„ 2 

/ 9-S^-dy [ y 2 h(y)dy, 



o Kv) 



o 



where we used Schwarz's inequality. The second integral 
on the right side is, integrating by parts, 



y 2 h(y)dy = / y 2 {~g'{y))dy = 2 1 yg{y)dy, 
o Jo Jo 



V. CONCLUSION 

We have studied the speed of travelling fronts of hyper- 
bolic reaction diffusion equations (^) for reaction terms 
for which the parabolic reaction diffusion equation (|l|) 
admits fronts. A variational principle for the speed of 
the fronts has been obtained for a wide class of systems. 
Explicit upper and lower bounds, obtained from the vari- 
ational expression allow one to characterize sytems for 



where we used the fact that g diverges at y = slower 
than (l/y 2 ) (this is equivalent to m/c > 1/2) and that 
g(l) = which can be seen form Eq.(l^). Replacing this 
in the inequality above we have the desired result, 



f yg(y)dy < 2 f g' 
Jo Jo 



(y)/h(y)dy. 
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lower bounds 


. f(u) = u(1-u)(1 +au) 

i < i i i < 


cKu) = 2-f'(u) 
i 1 i 1 i 1 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 

a 



FIG. 1: Bounds on the speed of fronts for an example 
with time delayed diffusion. The range of validity of linear 
marginal stability is reduced compared to the parabolic dif- 
fusion equation with the same reaction term. 



